bisection method with matlab code